Gauge theory approach to glass transitions 
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This theory combines a thermodynamic approach with a dynamic one in order to describe glass 
transition. Glass transition is regarded as an inaccessible second order phase transition, which is 
interrupted because of premature critical slowing down, caused by the system's frustration. The 
frustration-induced vortices are present in the structure besides thermoactivated vortices, and pre- 
vent the development of the order parameter fluctuations, that leads to the critical slowing down 
the system kinetics at some temperature above the phase transition point. 

PACS numbers: 



I. INTRODUCTION 

The formulation of the microscopic glass transition 
theory has remained one of the most intriguing but still 
unresolved problems of condensed matter physics for 
a long time. Many systems which manifest this phe- 
nomenon regardless of their nature enable us to con- 
clude that this phenomenon does not depend on any 
microscopic details, but it is determined by the symme- 
try properties of the systems, as in the case of phase 
transitions. On the other hand, this universality allows 
us to conclude that the glass transition is only the re- 
sult of general dynamic properties of condensed matter. 
Therefore, the question "Is the glass transition a phase 
transition, or is it the dynamic effect related to the lim- 
ited diffusion dynamics?" remains actual. One can con- 
cretize this problem by dividing it into some important 
questions that guide our presentation of glass-transition 
theories, which, for example, was formulated in [pi. In 
the present work we try to conciliate these positions, and 
answer some of the questions. 

We assume, that there are two key conditions for the 
implementation of glass transition. First of all, we be- 
lieve, that the glass transition has common nature with 
the second order phase transition, which is spontaneous 
breaking continuous symmetry and starting an ordering 
process in the system structure. On the other hand, 
in contrast to the phase transition this ordering process 
stops because of frustration, which arises in this process. 
The presence of frustration is the second key condition for 
glass transition. The frustration gives rise to production 
of vortexes (spin vortexes in the spin systems, or discli- 
nations in the undercooled liquids [ |, [| ) in the structure, 
which prevent the growth of the ordered regions. 

The nonequilibrium dynamics of the slowing vortex 
system, which takes into account the interaction between 
these vortexes, can be described in terms of the gauge 
field theory. Therefore the old, and almost forgotten idea 
of the gauge field description of glass transition underlies 
in this work j| [|, raj. Note, that this approach is very 
close ideologically with the "frustrated-limited domain 
theory" ||6|, [7), but the former seems to be more conve- 
nient versus the latter one, because of the absence of 
long-range interaction, which is deleted from the theory 



at the expense of the gauge field introduction. 

An important point of this theory is the presence of 
the vortexes in the structure. They can exist without 
frustration, but then their concentration tends to zero at 
the phase transition point, in which the relaxation time 
diverges. We take into account, that additional, induced 
by the frustration, vortices are present in the structure 
besides the thermoactivated vortices. Therefore, close to 
the glass transition there is a finite nonzero density of 
vortices. We believe that above the glass transition tem- 
perature the vortices system is in equilibrium. This fact 
makes it possible to average over the gauge field sources, 
that results in the temperature displacement of the crit- 
ical point in the gauge field subsystem. Then we use 
the renormalization group methods and critical dynam- 
ics methods for the description of the static and dynamic 
properties of the system at the new critical point. 

In terms of this approach the glass transition corre- 
sponds to the critical point in the gauge field subsys- 
tem. We examine both the static case and the dynamic 
case, and show that the theory's linear and non-linear 
susceptibilities agree with the experimental data. The 
dynamic features of this transition, such as the Vogel- 
Fulcher-Tamman relaxation time dependence on temper- 
ature, the characteristic temperature dependencies of the 
heat capacity, and the plateau in the time dependence of 
the order parameter correlation function, can be revealed 
by means of the functional methods of non-equilibrium 
dynamics. In conclusion we consider the relationship of 
the suggested theory with the mode-coupling theory ap- 
proach and the frustrated-limited domain theory, and 
discuss its physical interpretation. 



II. MODEL OF SUPERCOOLED LIQUID CLOSE 
TO GLASS TRANSITION 



One of the most elegant methods of description of the 
frustrated systems is given by the gauge theory g, ||| [| || ■ 
The static action of the system close to the second order 
phase transition in the general case has the form of 



S = P 



l( di s) 2 + U(s) 



dr. 



where s is some order parameter field, f3 = 1/TKb, Kb 
is the Boltzmann constant, and U(s) = /j 2 s 2 + vs 4 . The 
structure of vitrescent systems is inhomogeneous both 
in relative orientation of local ordering elements and in 
density. Congruent connections between the orientations 
are introduced by means of the gauge field, A^, when 
the ordinary derivative, 9iS, is replaced by the covariant 
derivative, Z^s, ||: 

t>iS a — > L>iS a = OiS a + g£abcAbiS C: 

where g is the Frank index (g rj 0.2 for a poly-tetrahedral 
packing), and e a b c is the Levi-Civita tensor. Besides, the 
system action should contain the gauge-invariant term 
~ (DiA a ) 2 . Therefore, one gets the gauge-symmetric 



action: 



S = /3 



^(D iS ) 2 + U(s) + ~F afiV F aiiV 



dr, (1) 



where 



As it was noted above, this model was suggested in [0, 0, 
II last century, but did not become widespread. 

In more detail the example of the three-dimensional 
Heisenberg spin-glass model, whose order parameter, s, 
is considered to be a local magnetization vector with the 
SO (3) continuous symmetry group, was considered in ||. 
With /x 2 > the action of this model is invariant under 
the SO (3) gauge transformations and (s) = 0. However, 
with /i 2 < the symmetry is explicitly broken, since the 
system can arbitrarily "choose" only one state from all 
equivalent states with the minimum energy U(s) poten- 
tial situated on the |s| = i[i/y2v sphere. We fix the vac- 
uum by means of fixing a point on the sphere. The system 
is no longer symmetrical with respect to the SO (3) gauge 
group, but it is invariant under the SO(2) group of the 
rotation around some chosen axis. 
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FIG. 1: Frustration induces topological defects (vortexes) in 
the system. These defects curve the space so that the local 
equilibrium configuration of spins contains the entire turn of 
their direction while going around this disclination along a 
closed contour. The availability of such vortexes in the system 
degenerates its collinear ground state. 

The expansion of the local magnetization field, s, 
near one of the vacuum states, for instance (s) = 
(0, 0, i(j,/\/2v), in small = s — ifi/^/2v deviations, 
and use of the gauge transformation properties allow to 



rewrite the action ([j]) in the form of the functional of two 
massive vector bosons, A Kfl (« = {1, 2}), with the mass 
Mq = ign/V2v, one massless vector boson, A^, and one 
scalar field, <j>: 



s = p 



(d^f + 2 M 2 2 



v<j? 



2 2 

g m 

4i> 



A A 



v A all A afi 



(2) 



dr. 



One of the important features of the above model (0) is 
the existence of soliton solutions, which correspond to the 
vortexes system (Fig.[l]). These vortexes are the movable 
sources of the gauge field, and their system above the 
glass transition temperature can be considered as a sys- 
tem which remains in thermal equilibrium [0. From (0) 
one can see that the gauge field correlation function di- 
verges in T c just as the scalar field correlation function, 
that corresponds to the sources concentration reaching 
zero at T+. 

The method of introduction of disorder into the the- 
ory is very important and plays a key part. In order 
for the system to possess glass properties we should in- 
ject disorder into this system, which gives rise to the 
frustration of its structure. In approach H the disor- 
der is associated directly with the quenched gauge field, 
and the A afJi field is frozen in an arbitrary configuration 
with some P(A afli ) distribution function. It is supposed 
that the gauge field describes the frustrations. But it is 
not quite right either, since the presence of the quenched 
gauge field in the system does not mean the presence of 
frustration yet. The frustration leads to multifold degen- 
eration of the structural states of the system. It produces 
the vortexes in the structure, which do not vanish in the 
low-temperature state. Therefore, one should introduce 
the frustration-induced vortexes into the model, in addi- 
tional to the thermally activated vortexes. 

The frustrated structure can be represented as vor- 
texes passing through the frustration planes [g, [| [Hj] 
(see Fig. I]). These vortexes are the sources of the A api 
gauge field, and should be injected in the model by means 
of some sources field, J api . Therefore, in contrast to j|, 
we believe that it is more correct to consider the source 
field, J aM , but not the gauge field, A afi , as the random 
field. In this case the A afJ _ field remains dynamic. For 
illustration the A afi field can be interpreted as a local 
relative rotation of neighboring local ordered domains, 
which corresponds to their local equilibrium. The do- 
main orientations can be movable, but their local equi- 
librium configuration around any frustration in the SV 
volume, bounded by the 5S sphere, should satisfy the 
following condition: 



2 f F aMV dS v 

ss 



J afJ ,dr ^ 0, 



(3) 



sv 



which follows from the principle of the least action. As 
a result one should add the term with the source field in 



the action: 



S = fi 



-(d^f + 2^ 



2j,2 



2 2 

9 V 

Av 



A A 



i . £ afiu F a^iv 



2 
v 4> + -^(t^Aa^Aa^ + J afJ .A all 



dr. 



(4) 

Let us consider the model (g|) in the fluctuation re- 
gion close to T c . Note, that this model was formulated 
for the low temperature state, T < T c , when the system 
symmetry is already broken. It is important, that in the 
high temperature state these reasonings are correct too. 
Mathematically it is formal extension of the theory to the 
high temperature region. One can understand the phys- 
ical meaning if the critical fluctuations are presented in 
the dual representation as moving vortices in the ordered 
phase. With approximation of T to T+ the density of the 
mobile vortexes decreases. The expression (0) binds the 
critical fluctuation scale with the inter-vortex distance 
and vortex density, so that the effective linear size of the 
critical fluctuations diverges proportionally to the gauge 
field correlation length. However, as it will be shown 
below, the frustration changes this picture a bit. 

At the temperatures above T c the structure is not a 
geometrical invariant, the vortexes (disclinations) can 
move, be born and be annihilated in it. But the frus- 
tration leads to the existence of a non-zero minimal vor- 
tex density. In terms of the field theory this means that 
any time in the system there is finite density, Iq, of the 
gauge field sources. Therefore the system can not be 
frozen in an ordered state. After freezing some vortices 
are quenched in the structure, but above this point the 
vortex subsystem is believed to be in thermal equilibrium 
H . Therefore, in the description of the state above the 
freezing point averaging over J afi leads to the redefinition 
of the partition function: 



Z = 



exp 



--! 



I l J 2 ail dv ) SJ a , 



®<j>®A afJ 



freezes in a disordered state. The T g can be regarded as 
the glass transition temperature. We analize this suppo- 
sition below. 



III. ANALYSIS OF THE SPIN-GLASS MODEL 
IN STATIC CASE 



As it was noted above, averaging over J^ leads to shift- 
ing the transition temperature from T c to T g . In this case 
the mass of the gauge field becomes M, and the propa- 
gators of the theory fields are written in the following 
form: 



Go(fc) = 



1 



r + M 



A^ = 



k 2 + M 2 



Let us choose the model parametrization so that T g Ks = 
1 for simplicity. Close to the glass transition M —> 0, 
and the gauge field correlator becomes the same as the 
propagator of a massless photon, 

A -^ 
"" ~ k 2 ■ 

It is known that the theory is renormalizable in this case. 
Therefore, one can investigate the critical properties of 
the model at T ^ T g + . 

The canonical dimensions of the values of our theory 
are given in the table: 



F 


k 





A 


M 2 


M 2 


.9 


V 


d k [F] 


1 




i rffc 


4-2 


4-2 









The renormalization procedure only refines these values, 
which results in the replacement of the canonical dimen- 
sions by the critical ones A[F) = d[F] + 0[F](e), where 
e = 4 — dk, and dk is the space dimension. 



where J . . . T)x is the continual integral. It results in the 
additional contribution to the A afl "mass", which is as 
follows: 



M 2 = M 2 -k ) = n 2 g 2 /Av-h. 



(5) 



Thus, the frustration leads to the renormalization of the 
gauge field mass. 

The renormalization of the gauge field mass affects the 
critical behavior of the system, since it shifts M 2 = 
singularity to the temperature range above the virtual 
phase transition point, T c . If we assume that [i 2 = 
aKsiT — T c ), where a is some constant, then from (H) 
we have the critical divergence of the A a field correlation 
radius at T g = T c + AIov/aKsg 2 - In g] it was shown 
that it results in the critical slowing-down of the fluctua- 
tions. Thus, the disorder-induced frustrations inhibit the 
growth of the cf> field correlation length, and the system 
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FIG. 2: The graph representation of the Green functions: 
a — Go, b 



A,j„, and the action terms. 



The renormalization procedure is carried out with 
the standard method prf. It is assumed that the 
fields <p an d A are slow- varying ones, such that the 
Fourier-transformed fields have only long-wave compo- 
nents: |fc| < fco- At the first step of RG transformations 
one integrates the partition function over the components 
of the fields in the limited wave band A/co < k < ko, 
where A <$^ 1 is the regularization parameter (momentum 
cutoff). It is expected that under certain conditions this 
action has a structure similar to the original one, in this 



case a model is multiplicatively renormalizable. One can 
check that the formulated model satisfies this criterion. 
As a result one gets an effective action S\ with renormal- 
ized parameters (Z v , Z M , Zm, Z g , and Z g i), which are 
named the "constants of renormalization" and depend on 
the cutoff A; at the second step, one makes the inverse 
scaling transformation of the fields (Za, Z$) and coor- 
dinates, which aims at restoring the original cutoff scale 
fco. Then the renormalized parameters have the following 
form: 



M 2 ^=Z M ,Z A A d = A- 2 Z M ,, 
A- E/2 Z n , 



gW = Z g Z A A 2d+1 

g 2(R) = Z 2 Z A A 3d = 

vW = Zv ZiA 3d 



(6) 



= Z 2 Z 2 A Z 2 A 3d+3z 

A- £ Z„. 



A~ e Z, 



tr' 



It is assumed that S is invariant with respect to the above 
scale transformations at the critical point. 

Close to M 2 = (T > T g ) the gauge field be- 
comes massless, but the scalar field remains massive with 
/j 2 = 4Iov/g 2 . Therefore, the contribution to the renor- 
malization is made only by the loops of the gauge field 
propagators. 

In Fig. y some graphs giving logarithmically divergent 
contributions to the renormalized theory are presented. 



where £ = ln(l/A) is the logarithmically divergent fac- 
tor. In the same way one can get other terms of the 
renormalized action: 



g W = e^ 2 Z g ~ e^ 2 

g 2(R) = e eZ Zg2 „ e ei 

V W = e*Z v ~ e e « 



2 MY 



8^ 



9 ~^ 
2 .9% 



(9) 



Hence, in the one-loop approximation the renormaliza- 
tion group of the model under study has the form: 



<91n(M 2 ; 
d\n(fi 2 ) 

31n(g 2 ) 
d\nv 



= 2 - 3.g 2 /87T 2 , 

„ M 2 g 2 „ 

8 M V ' 



(10) 



£ - .9 2 /4tt 2 



c>£ 



e-g 4 /Sv-K 2 . 



a 

X'"'X 



From the condition of the stable point existence, 
<91n( ff 2 )/d£ = 0, dhx(v)/d£ = 0, we get 3 2 = 4tt 2 £, and 
v = <? 2 /2. In addition, one can see that 
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FIG. 3: The graphical representation of one-loop contribu- 
tions to the terms of the action. 

Let us explain in detail the renormalization of M 2 as 
an example. We will limit ourselves to using the one loop 
approximation, that is quite enough for the demonstra- 
tion of all the features of the theory. In this case the 
renormalization constant of M 2 has the form 

k Q 

6M 2 g 2 



Zm^S^u ~ M 5^ v 



4(27r) rffc 



M 2 S 



3M 2 g 2 S 



uvUkir 



, , A M (fe)A A „(fe)dk 

dfc/2 



/ ' ' ' 



(2ir) dk r(l + d fc /2) 



fc dfc 



-dfc. 



Afcn 



One can see that the integral in this expression introduces 
a logarithmically divergent contribution to M 2 renormal- 
ization if the momentum dimension is dk — 4. In this case 
we get following expression for the renormalized value of 
M 2 : 



M 2(R) = e 2f z 



M 2 



■M 



M 2 -3^' 



(8) 



M 1 = 



aKsg 



(T - T g ) « e 



2S 



(11) 



Thus, in T g the gauge field becomes massless, whereas 
the local magnetization field still has got the "mass" , /i. 
The analysis of the model close to this point allows us 
to assert that T g is nothing else but the glass transition 
temperature. First of all, it is evident from (0) that the 
linear susceptibility, \l = d{<j>)/dh ~ /i~ 2 = g 2 /AIqv (h 
is an external source of the field 0), is finite at T = T g . 
Similar unsophisticated estimation gives the correlation 
length, r cor ~ y/ g 2 /AIqv, is finite too. One can see that 
in the absence of the frustration, Iq = 0, the linear sus- 
ceptibility and the correlation length diverge in T g . This 
is natural, because in this case T g = T c and the glass 
transition becomes the second order phase transition. 

On the other hand, nonlinear susceptibility, \n — 
d 3 (cj))/dh 3 = (0 4 ) p2[ , diverges near T g because of the 
divergence of the A-field loop contribution (Fig. Eh: 



XN 



^ ~{T-T g 



(12) 



where 7 sa e/(2 — 3e/2). The estimation of this exponent 
in |9| was wrong because of incorrect summation of the 
graph series. The more accurate estimation gives 7 sw 2 
when e = 1. 

Using this theory allows us to determine the tempera- 



ture dependence of the system heat capacity at T —} T + : 



d JL^-K B \nZ- 



p dT 
Vag 2 

+—r- 



K B TdZ 
Z dT 



-KnlnZ 



^-i)[ir-T,>-iw.i<,v-o" 



rn| A_i.) rr t, )■,>,,:,.-„] 



+3V z v 



d 



1 



dT T 



(13) 



[(#)r=o(#)r=o] 



2 / r) 1 

+y2 ^4~(aT^^' [(.l;,-J;y).-o(-i;,---l;;).- o] 



Since the matter keeps disordered structure at freezing, 
we will believe that Z weakly depends on the temperature 
at T g . Let us assume that the first term is some constant. 
One can present the correlation functions in the form of 



\AijA. 



%J- n *3lr 



= H -l e -k|M 



(ipip), 



= | r |-l e -Hr 



(14) 



where M ex (T-TJ 1 / 2 , and m ex (T-T c ) x l 2 . Therefore, 
close to T g the terms with J^{AijA i:j ) oc (T - Tg)' 1 / 2 
are dominant, since they diverge at T g , as opposed to 
the terms which contain only (tptp). As a result c p ex 
(T-T s )-V2. 

The above results show that T g is a critical point of the 
model, but they do not allow to conclude that it is the 
glass transition that occurs at T g . The correlation length 
of the ^4-field diverges in T g , but the physical meaning 
of this correlation function is not clear. Besides, the ob- 
tained temperature dependence of the heat capacity near 
T g diverges and is not characteristic for glass transition. 
One can show that this is a result of only the static char- 
acter of the theory. Below the analysis of the critical 
dynamics of this system close to T g allows to make sure 
in it, and to show that T g is the glass transition temper- 
ature. 
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FIG. 4: The qualitative presentation of the temperature func- 
tion of the nonlinear susceptibility, xn oc 1/(T — T fl ) 7 , and 
linear susceptibility, xl oc 1/(T — T c ). The xl curve is not 
described by this theory at T < T g , it is finite in this region. 
The qualitative picture of the temperature dependence of the 
heat capacity near T g , which is obtained from the theory. 



with 



S* 



$(t,T)G- l (t-1?,T-T>)$V,T>) 



drdr'dtdt' 



+A a ^{t, p)A-^(* - 1', r - r')A av (t', r') 

[ge a bc(d^A au )A hl A cu + ge abc (d^A a u)A bll A cll 

+ g£ abc\0 \iA av ) A b ^A cv + g S abc^ aij Ab^A cv Ai^Aj v 



IV. SPIN-GLASS MODEL IN DYNAMIC CASE 

One can examine the non-equilibrium dynamics of the 
system close to T g by means of functional methods of 
non-equilibrium dynamics |L3|| near the critical point. It 
leads to the representation of the partition function of 
the system in the form of 



Z = / exp(-£*)2)0D4». 



/'•• 



(15) 



+g 2 A a „ A a ^ 2 + g 2 (A at ,) 2 cj)(j) + v4 # 3 ] drdt, 

( 16 ) 
where = {</>, </>}, and A a{1 _ = {A api , A api } are vec- 
tors, the components of which are named "quantum" 
and "classical", respectively, in the Keldysh representa- 
tion |13[| , G _1 and A~^ are matrices, inverse to the Green 
functions matrices having the following form: 



G = I °„ ° L A, 



G$ 



A K A A \ 

a r o r [ ' 



XXX 



V 



e . f , 

■ 
A ifi 



■ 



FIG. 5: The graph representation of (ha) action terms, 
where 



k + fi ± ir^w 



G*(fc, W ) 



2r, 



(fc 2 + Ai 2 ) 2 + r 2 u, 2 ' 



(18) 



r^ is the kinetic coefficient of the local magnetization. In 
the case of M -> 



A^(fc,u;) 



J liu 



k ± ir^w ' 



71 i ti2 2 : 
K + 1 > 



A^(fe, w) 



(19) 



where T^ is the kinetic coefficient of the gauge field. 

a b c 



where 4[-F] is the frequency dimension [[III Q. The 
canonical dimensions of the values of our theory are given 
in the table: 



F 


k 


LO 


4> 


^ 


A 


A 


d k [F] 


1 





i 4 

-!-T 


i 4 


i 4 


1 4 


d u [F] 





1 


-1 





-1 





d[F] 


1 


z = 2 


o 4 


i d k 


o 4 


i 4 



F 


.9 


u 


T A 


r* 


M 2 


M 2 


4[F] 








4-2 


d k -2 


4-2 


4-2 


4[F] 








-1 


-1 








d[F] 








4-4 


4-4 


4-2 


4-2 



The renormalization procedure only refines these values, 
that leads to the replacement of the canonical dimensions 
by the critical ones A[F] = d[F] + 0[F](e). 

The renormalization procedure is carried out with 
the standard method. It is assumed that the fields 
<f>, (f>, A, and A are slow-varying ones, such that the 
Fourier-transformed fields have only long-wave compo- 
nents: \k\ < k ; uj < cj . At the first step of RG 
transformations one integrates the partition function over 
the components of the fields in the limited wave band 
Afco < k < kg, A z ujq < w < luq. The renormalized pa- 
rameters have the following form: 



FIG. 6: The graph representation of the Green functions: 
a - G$, b - Gt c - G?, d - A* , e - A^, f - A* . 



V. CRITICAL DYNAMICS OF THE MODEL 

The critical behavior of the system close to the T g 
can be considered within the critical dynamics tech- 
nique 1^, Q. The critical dynamics rests on the hy- 
pothesis of dynamic scaling, according to which the ac- 
tion should be invariant under scale transformations, 
which conformly expand the space and time coordinates 
(u ~ k d ^). In this case the summarized dimension, 
d = dk + 4 (4 = z is the dynamic exponent), has 
the same role as the conventional (momentum) dimen- 
sion, 4i in the static case. The canonical dimensions of 
the fields and the model parameters are determined from 
the condition of dimensionless action. The corresponding 
summarized canonical dimensions, d[F], of any values, F, 
depend as 



/?(«) = VZ^A d +*= A- 2 V, 
M 2 ^ = Z M 2Z A Z A A d+z = A- 2 Z M 2, 
rW = Z rA Z A ZzA*+' = A e Z rA , 
qW = Z ZzZ 2 A A 2d+2z+1 = A~ e / 2 Z, 



(20) 



g 2(B) = z 2 g Z A Z\ A 3d+3z 



(/■■ 



A~ e Z, 



v 



2 ^ = Z 2 Z A Z A Z 2 A 3d+3z - 
{R) = Z„ZzZlA M+3z = A" 



9^ 

: A~°Z g 

l z.„. 



d[F] = d k [F] + z-d ul [F} 1 



As well as in the static case, close to M 2 — (T ~ T g ) 
the gauge field becomes massless, but the scalar field re- 
mains the massive with /i 2 = Alov/g 2 . Therefore, the 
contribution to the renormalization is made only by the 
loops of the gauge field propagators. According to the 
separation of massive field theorem |15| the Feynman di- 
agrams, containing the propagators of the field, the mass 
of which is appreciably larger than the external momen- 
tum, are inversely proportional to the degree of this mass, 
and make a finite contribution to the renormalization. In 
Fig. u\ some graphs giving logarithmically divergent con- 
tributions to the renormalized theory are presented. 



In the one-loop approximation the renormalization 
constant of M 2 has the form 

Zm 2 S^ ^ M 2 5 tll/ 



(2tt 



A/,',, 



fe0 _J^ -^fc*-l 



Af 2 <W 



12MV /" r(i + 4/2) 



(27r)^ +i J (k 4 + co 2 )(k 2 + iuj) 

Ak 



-dkduj 



where s = 4 — d& . 

One can see that the considered theory is renormal- 
izable, since the renormalization procedure leads to the 
correction of only existing terms of the Lagrangian, and 
all renormalized parameters logarithmically diverge when 
dk = 4, which thereby is the critical dimension. 

The interaction of the local magnetization fluctuations 
per gauge field (Fig.H) plays the key part in the con- 
sidered theory. This becomes clear when we consider 
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"\M 2 n 2 A d,-rr dk / 2 f k dk ~ 1 
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Afe Q 

The integral in this expression introduces a logarithmi- 
cally divergent contribution to M 2 . In this case the ex- 
pression for the renormalized value of the M 2 coincides 
with (§1): 
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FIG. 



in detail the contribution to the renormalization of the 
r^ node of the a diagram, which is shown in Fig.ffl. 
This term is also interesting because the T^ node is pro- 
portional to the relaxation time of the local magneti- 
zation field, and determines the kinetics of the glass- 
transition: the loop of the light field, A a : which is given 
in Fig.g, makes the logarithmically divergent contribu- 
tion ~ ln(l/A)<5(u;). This term determines the divergence 
contribution of the a graph (Fig. 0) in which the mas- 
sive field loop leads only to multiplying the logarithm by 
4I v/g 2 factor: 
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FIG. 7: The graph representation of the one-loop contribu- 
tions to the terms of the action. 

In the same way one can get other terms of the renor- 
malized action: 
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(22) 



where t is the time of the observation of the system. 
One can see that for I —> 0, or with a short observa- 
tion time, T r f,g 2 /4lQV 3> t , this contribution becomes 
negligibly small, and the theory becomes nonrenormal- 
izable. This relates to the symmetry properties of the 
Yang-Mills model. In this case the fluctuation dissipa- 
tion theorem (FDT) is always broken because of the free 
energy transfer between the local magnetization modes 
and massless gauge field modes (Goldstone modes). How- 
ever, there is also some problem in the presence of the 
quenched disorder: this node becomes nonlocal in time. 
This formally makes the renormalization group equa- 
tions nonautonomous for t « T,pg 2 /AIqv, which violates 
the correctness of the renormalization procedure. The 
nonautonomous renormalization group has recently been 
considered in physical problems rather often (see, e.g., 
fll6| , |l7f). In order to avoid this problem we note that 
the renormalization group equations are autonomous in 
two limit cases: t <C T^g 2 /AI v and t ~^> r '^g 2 / 'AIqv , 
the physical solution is the result of matching these two 
limit cases. 

At large time scales, t 3> V ^g 2 / AIqv : the contribu- 
tion to renormalization of T^ is logarithmically divergent. 
Hence, in the one-loop approximation the renormaliza- 



tion group of the model under study has the form: 
dln(I» 
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(23) 
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As well as in the static case, from the condition of the 
stable point existence, <91n(g 2 )/<9£ = 0, <91n(u)/9£ = 0, 
we obtain g 2 = 4ir 2 s, and v — g 2 /2. 

At small time scales, t <C T^g 2 /AI v, the contribution 
to the r^ renormalization is negligibly small. In this case 
the RG-equations form similarly to (E3J) except the first 
equation, which has the form: 






= 0. 



(24) 



In order to match the solutions found in ( |23| ) and (E4) 
it is necessary to analyze the renormalization in the t ~ 
T r j > g 2 /AI v region, and get a matching function, $(£), as 
it was done in fi \\M : 






5 4 $(0/7T 2 



where £ ~ 1. The matching function, which is 1 in the 
case of large time scales, £ ^> 1, and in the case of small 
time scales, £ <C 1, can be chosen in the form: 

$(0 = 1 - K z = 1 - exp(-*0, 
where z ~ 2 is the dynamical index O, ttm . As a result 
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(25) 



Hence, critical slowing down of all relaxation pro- 
cesses does occur at T g , and follows the Vogel-Fulcher- 
Tammann relation, which in our case was derived from 
the microscopic reasons by means of the Keldysh tech- 
nique and critical dynamics method. 

Note that if frustrations are absent, Iq — > 0, the freez- 
ing temperature coincides with the phase transition tem- 
perature. Then the diagrams with the loops of </> and 
4> fields become divergent, and the system experiences 
the paramagnetic-ferromagnetic phase transition, which 
is described within the standard critical dynamics fll9| . 

In the one-loop approximation the (/>-field correlation 
function can be represented as the sum of the unper- 



turbed and cooperative parts: 



2I og \HL/a ) J G K {t , )e -u a «\ t -t>\/r^ d ^ (26) 
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where L is the size of the system, and a® is the inter- 
atomic distance. This function has the form which is 
characteristic for glass systems (Fig. ft). 




ln(t) 



FIG. 9: The dependence (</></>) t on ln(t): a) is the contribution 
of {^><j>)t, which is the Debye-relaxation; b) is the contribution 
of the second term which is given by the cooperative effects; 
c) is the sum of the first and second terms. 



The approximation of the linear and nonlinear sus- 
ceptibilities in the dynamic theory coincides with that 
the static one. However, the dynamic theory allows us 
to improve the temperature dependence of the system 
heat capacity at T — > TT. Carrying out the analy- 
sis, it is convenient to consider the correlation functions 
in the (t, ^-representation, and take into account that 
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0. As a result we obtain: 
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where the /,(...) = J_ {■ ■ ■ )dk notation is introduced. 
Since the matter keeps disordered structure at freezing, 



we will believe that Z weakly depends on the temperature 
at T g . Therefore, let us assume that the first term is some 



constant. From (J18|) and (19) it is not difficult to get: 
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where £/c(to) = k 2 + to 2 , whence 



-\t\s k (M)/r A (T) 



(27) 



(28) 



1 



t=o — 



VT V (TY J k 
Consequently, we obtain: 
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Since the system is close to T g , then the last term gives 
the most considerable contribution in the heat capacity. 
One can obtain the qualitative form of the temperature 
dependence of the heat capacity close to the glass tran- 
sition: 



c p (T) ex 



d_ 

df 



-CT g /(T-T g ) 



-CT g /(T-T g ) 



(T-T g 



-CT g , 



where C — 2vg 2 /air 2 . The graph of this expression is 
shown in Fig.Q. In this picture one can see the sharp 
growth of the heat capacity near the glass transition 
which then drops at T gi which is characteristic for vitres- 
cent systems. 

According to the suggested theory one can give a sim- 
ple qualitative explanation of this dependence: when 
the system temperature approaches at T c , the scale of 
the thermal fluctuations increases, and consequently, the 
heat capacity increases too, as in the case of the second 
order phase transition. However, critical slowing down 
inhibits the growth of these fluctuations. The maximal 
size of the them is determined by the frustration scales. 
As a result the heat capacity does not infinitely diverge 
at T g , but gets a finite maximum at T m ax > T g , and 
then falls down (Fig.[|). 



VI. DERIVATION OF THE MODE-COUPLING 
THEORY EQUATION 

As the suggested gauge theory is based on the micro- 
scopic conceptions and has a general character, so it is 



natural to trace its relation with other known theoreti- 
cal approaches to the description of glass transition. The 
most applicable of them is the Mode-Coupling Theory, 
which, in the Zwanzig-Mori representation has the form 
of the motion equation for the pair correlation function. 
The key problem of this theory is the determination of 
the memory function, which is the kernel in the integral- 
differential equation. 

In order to derive the motion equation for the model 
structural factor, S{t) = Gjf(£), in the presented theory, 
one can use the Dyson equations for G A and G R : 



C A = C A ^ 



+ G A{R) ®D A ^®G A ^ 1 (29) 



where ® denotes the convolution of two functions on t, 
and 



D A(R \t) = 6{±t) 



g yin(£/ao) r *4i n vt/r^ 



2T rh TT Z 



(30) 



is the self-energy part, which has the physical meaning of 
the memory function in the mode coupling theory. After 
the Gq 1 — T,pdt — V 2 + fi 2 operator action on these func- 
tions one can subtract the first equation from the second 



G^iG* - G A ) = D R ®G R -D A ®G A 
t 

D R (\t'\)G R (\t~t'\)dt' 

o 



D R {\t'\)d t G K {t-t')dt'. (31) 



Then, using the FDT, one can obtain: 



d 2 G K (t)+^ 2 d t G K (t)+T / D R (\t'\)d t G K (t-t')dt' =0. 

o 

In order to present this expression in the typical form, 
one should add the static part to the Green function and 
memory function, D(t) — > D{t) + Vi 2 jT . $7 is the micro- 
scopic frequency, obtained from the experimental static 
structural factor. Then we have: 

d 2 G K (t) + n 2 G K (t) + ti 2 d t G K (t) 

i 

+ T f D R {\t'\)d t G K (t - t')dt' = 0. 
o 

This equation being exact, the crux of the mode-coupling 
approach consists in formulating an approximate expres- 
sion for D R {\t'\). In the considered theory this function 
can be approximated by the sum of the diagrams which 
give contribution to the self-energy part of the Green 
function. 
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VII. RELATION WITH THE 
FRUSTRATED-LIMITED DOMAIN THEORY 

The suggested theory is related very closely both with 
the Stillinger's "tear and repair" mechanism for relax- 
ation |2(| and the "frustrated-limited domain theory" of 
Kivelson and Tarjus []6|. In these theories the frustration 
is described as the source of the strain free energy that op- 
poses the spatial extension of the locally preferred struc- 
ture and grows with the system size. It results in break- 
ing up the liquid structure into domains, whose sizes are 
limited by frustration with decreasing temperature. 

One can switch over from the model (f|) to the 
frustrated-limited domain theory. Let us not average over 
J-field, as it was done above, and suppose that the discli- 
nations are arbitrarily located in the structure. Instead 
of this let us get rid of the gauge field and carry out the 
integration of ([|) over A a . In order to simplify the con- 
sideration one can make a rather strong approximation, 
that the system relaxation depends only on the kinetics of 
the vortex system, and one can neglect the non-quadratic 
terms in (Q). Then the functional integration over A®, 



Z A = / ®A a u exp 



\a%x£ai + j%a% ) ,/k 



leads to 



Za oc exp 



, J a u k~ 2 J a u dk 



exp 



i|j M Q (r)|r-r'rV M V)drdr' 



As a result, the system action can be reduced to the 
following form: 

S = S,p + Sj, 

where S^ is the scalar field part of the action, and 



sj^c^j^j;^ 



c, 



E 



JSWufa) 



where the disclination element J?i{ T i) occupying the site 
i at r, position, is the contribution of the geometrical 
frustration M . The first term is the disclination core con- 
tribution to the action, and the second term is the con- 
tribution of the vortexes (disclinations) interaction [El] . 
According to the Kivelson-Tarjus theory the competition 
of these terms determines the properties of glass-forming 
liquids close to the glass transition. 



VIII. CONCLUSIONS 

According to the above theory one can offer the fol- 
lowing physical picture of the glass transition processes: 
when the temperature approximates T g , the correlation 



length of the gauge field diverges, while the correlation 
length of the order parameter is finite. The infinite corre- 
lation of the gauge field means that the relative rotations 
arc correlated at an infinitely large distance. Then the 
growth of the order parameter field fluctuations becomes 
impossible, since the spins can not already turn indepen- 
dently of each other. 

The statement about that the correlation length of the 
gauge field grows faster than the correlation length of 
the order parameter can be explained in the following 
way: without the frustrations the correlation lengths of 
both the gauge field and the order parameter field grow 
equally, because they are components of the same field. 
In the dual representation the system contains thermally 
activated vortexes, whose concentration tends to zero 
when T — > T+. The priority growth rate of the correla- 
tion length of the gauge field becomes possible when the 
system contains frustration. The frustration induces vor- 
texes additional to the thermally activated vortexes, as a 
result in the equilibrium state the vortex concentration 
is fixed non-zero. The reason for this is that the fluctua- 
tions of the gauge field are sure to develop around these 
sources. In the case of a nonzero concentration of these 
static sources the fluctuations around them can interflow, 
which leads to the faster growth of the effective correla- 
tion length of the gauge field. With the given source con- 
centration and the appropriate correlation length of the 
order parameter (or the gauge field without frustrations) 
it leads to the formation of the percolation cluster which 
is associated with these sources. Thus, the effective cor- 
relation length of the gauge field becomes infinitely large, 
and the relative rotations at infinitely remote points be- 
come correlated and the relaxation time of the system 
infinitely grows up, but the order parameter correlation 
length remains to be relatively small. Therefore, the sys- 
tem freezes in a disordered state. 

In conclusion, note that the above theoretical approach 
describes all general properties of glass transition very 
well: the expression for the temperature dependence of 
the heat capacity near the glass transition is in good 
qualitative agreement with the experimental data. The 
derived temperature dependences of the linear and non- 
linear susceptibilities behave quite properly at T g , and 
the critical exponent of the non-linear susceptibilities is 
7 ~ 2, that is a good estimation. To crown it all, this 
theory reproduces the characteristic form of the {4><t>)t 
correlation function dependence on time, and predicts 
the Vogel-Fulcher-Tammann law (E5h for the tempera- 
ture dependence of the relaxation time. The theory al- 
lows to answer some important questions that guide our 
presentation of the glass-transition theories JjJ: first of 
all, according to the suggested theory, the time of relax- 
ation increases when T — > T+ with non-Arrhenius law 
( p5| ) because of the features of the non-equilibrium dy- 
namics of the frustrated system in the critical fluctua- 
tion region. The broad relaxation spectra are the result 
of the dynamic scaling in the fluctuation region, which 
results in the hierarchy of relaxation times. It is very 
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difficult to say something about the relation between the 
kinetics and thermodynamics discussing the glass transi- 
tion within the presented theory. One can only assume 
that T c can be considered as the Kauzmann temperature, 
T c = Tk, because both of them are virtual parameters 
with similar meaning. In T c the non-frustrated system 
tends to the hypothetical ordered state with the mini- 
mum of the configuration entropy, but this state, as well 
as the Kauzmann state in T = Tk, is not physical. So 
the presented theory is based on the suggestion, that in 
the system there is an inaccessible (virtual) transition, 
T c , but this is not an "ideal glass transition", since this 
inaccessible phase is assumed to be ordered. The collec- 
tive motion in the theory can be represented as a vortexes 
motion, therefore the glass transition, when this motion 
becomes correlated and critically slows, is a collective 
phenomenon. The term "order parameter" has somewhat 
different meaning in comparison with the phase transi- 



tion theory. Here the "order parameter" is the material 
field, which corresponds to the local ordering with the 
structure of the hypothetical ordered state as the local 
magnetization in spin systems. It is the order parameter 
of the inaccessible (virtual) phase transition, T c . Accord- 
ing to the suggested theory the dynamic correlation is di- 
rectly associated with the gauge field correlation, because 
it does determine the divergence of the relaxation time of 
the order parameter. This theory gives a good geometri- 
cal representation of the glass transition as freezing the 
vortices network. It was developed for the spin systems, 
but it can also be reformulated for the description of the 
supercooled liquid-glass transition. In both cases it has 
the same topological nature, therefore, one can expect 
qualitatively identical results. In this sense the theory 
has a general character. 
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